Quantum and Classical Glass Transitions in LiHoa;Yi_a;F4 
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When performed in the proper low field, low frequency limits, measurements of the dynamics and 
the nonlinear susceptibility in the model Ising magnet in transverse field, LiHoxYi_xF4, prove the 
existence of a spin glass transition for x — 0.167 and 0.198. The classical behavior tracks for the two 
concentrations, but the behavior in the quantum regime at large transverse fields differs because of 
the competing effects of quantum entanglement and random fields. 
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Research on spin glasses Q has led not only to deep 
insights into disordered materials and the glassy state, 
but has generated novel approaches to problems rang- 
ing from computer architecture through protein folding 
to economics. The rugged free energy landscape charac- 
teristic of such systems defies usual equilibrium analyses, 
with pronounced non-linear responses and history depen- 
dence. At low temperatures, and in cases where barriers 
to relaxation are tall and narrow, quantum mechanics can 
enhance the ability to traverse the free energy surface Q . 
The LiHoa^Yi_a;F4 family of materials represents the sim- 
plest quantum spin model, the Ising magnet in transverse 
field, and it has been an especially useful system to probe 
the interplay of disorder, glassiness, random magnetic 
fields and quantum entanglement [1, 0, [l, 0, 0, H, 0, [l^ • 
The parent compound, LiHoF4, is a dipole-coupled Ising 
ferromagnet with Curie temperature Tc = 1.53 K. Ap- 
plying a magnetic field Ht transverse to the Ising axis 
introduces quantum mixing of classical spin-up and spin- 
down eigenstates which are split by an energy F ~ H^^ 
or equivalently, tunes the tunneling probabilities for walls 
between patches of ordered spins [11]. Hence, quantum 
fluctuations controllable by an external field can drive 
the classical order-disorder transition to zero tempera- 
ture, resulting in a much studied ferromagnetic quantum 
critical point pj, UJ, [jj • 

The nature of the ground state can be tuned by par- 
tially substituting non-magnetic Y for the magnetic Ho 
P^. Dilution enhances the effects of the frustration in- 
herent in the dipolar interaction, with the ferromagnet 
giving way to a spin glass at x ~ 0.2 (Fig.[T]). The trans- 
verse component of the dipolar coupling introduces other 
phenomena. At large x, the ground state becomes a dis- 
ordered ferromagnet, with the low temperature dynam- 
ics dominated by domain wall tunneling [ll|. The high 
T and low Hf behavior reveal both the effects of Griffiths 
singularities [16] and the internal random fields due to the 
application of a uniform Ht to a disordered Ising mag- 
net 0, H, ^ ^ 0-1 5 internal transverse fields, 
F^, induce quantum entanglement that prevents the sys- 
tem from freezing and stabilizes a spin liquid "antiglass" 
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FIG. 1: (color online) Energy levels and phase diagram for 
LiHoccYi_ccF4 (a) Schematic of energy levels for antiferromag- 
netically coupled spins in a (top) random field and (bottom) 
uniform transverse field F, and then a random field h. (b) 
Magnetic phases in the x — T plane. Arrow denotes spin liquid 
"anti-glass" phase, (c) Magnetic phases in the x — Ht plane. 
Open circle shows peak in x{Ht) for the antiglass [l8|. (d) 
Spin glass/paramagnet phase boundaries for x = 0.167 and 
0.198 cross. 



phase down to very low temperatures [l, 0, [l^ . In the 
intermediate range, where the spin glass phase is sta- 
ble, the tendencies towards ferromagnet ism and random 
field effects found at high x compete with the massive 
quantum entanglement of the antiglass. This competi- 
tion might be expected to lead to very different statics 
and dynamics for small changes in x. Fig. [1^ illustrates 
the essential physics. In the classical (F^ = 0) case, a 
small "random" field of strength h at site 1 and at site 
2 will produce a splitting of order h between the degen- 
erate I ti)^ I it) classical ground states. On the other 
hand, Ht will yield a single non-degenerate ground 
state I II) + I II), on which the only effect of a small h 
will be an energy change ~ h'^. 

The first work on the spin glass state in 
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LiHoo.i67Yo.833F4 revealed a well-defined transition 
from the paramagnet to the glass, evidenced by both a 
sharp divergence in the nonlinear susceptibility, xs^ and 
a dynamical signature in the dissipative component of 
the linear susceptibility, ^' [sl, (3] . These measurements 
relied on Rt to speed up the system dynamics to the 
point where the / ^ limit could be probed directly, 
and led to considerable theoretical work 0, H [ii, 
some with good qualitative agreement with our experi- 
ments [20]. Within the last three years, the concept that 
randomly placed classical dipoles [2l] should undergo 
a spin glass transition has itself been questioned on 
account of numerical work on small cubic lattices [22]. 
More recent experiments [23| employing a /iSQUID 
magnetometer did not reveal a divergence in xa, lead- 
ing the authors to a similar conclusion, namely that 
LiHoo.i67Yo.833F4 is not a spin glass. Unfortunately, 
the authors of [23| used large longitudinal fields and 
fast sweep rates, probing the system very far from 
equilibrium, and thus obscured the meaning of their 
data. In this Letter, we show explicitly that when data 
are acquired in the proper small longitudinal field, low 
frequency limit, clear evidence is seen for a spin glass 
transition. Moreover, we report the discovery that a 
minor change in x from 0.167 to 0.198 - approaching the 
multicritical point where spin glass and ferromagnetic 
phases coexist (Fig. [TJd) - results in dramatic changes in 
the quantum (i^^-dependent) behavior. 

We performed ac susceptibility measurements from 
1 to 10^ Hz on single crystal needles of LiHoa;Yi_ajF4 
mounted on the cold finger of a dilution refrigerator. The 
magnitude of xi was consistent with the Curie- Weiss law 
for Ho^+ ions at high T. The susceptibility was corrected 
using the demagnetization factor of rods of the same 
aspect ratio. Ho concentrations x were determined to 
zbO.OOl by a differential weighing technique. Static trans- 
verse magnetic fields up to 80 kOe and longitudinal fields, 
/i^, up to 300 Oe were supplied by a superconducting 
solenoid and Helmholtz coils, respectively. The ac exci- 
tation amplitude was restricted to less than A = 0.02 Oe 
to ensure linear response and to control heating. For 
measurements of the nonlinear susceptibility, the dc lon- 
gitudinal field was swept at 0.04 Oe/s so that this rate 
was smaller than the effective sweep rate at /= 1.5 Hz 
of 27rA/ = 0.2 Oe/s. 

We plot in Fig. [T]i the T-Rt phase diagrams for both 
the X = 0.167 and 0.198 spin glasses. The transition is de- 
fined by the emergence of a flat spectral response at low 
/ in Xi , corresponding to 1// noise in the magnetization. 
This dynamically-determined phase boundary coincides 
with that derived from the maxima of X3(/ ~^ 0)- The 
classical spin glass transition Tg{Ht = 0) increases with 
increasing x, but lies below the mean-field ferromagnetic 
Tc{x) = xTc{x = 1). However, once the transverse mag- 
netic field is turned on and the relative importance of 
quantum entanglement and random field effects becomes 
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FIG. 2: (color online) Scaled susceptibility for x = 0.167, 
Ht = 2.1 kOe (a) and x = 0.198, T = 0.25 K (b) as a function 
of longitudinal bias field added to sinusoidal 1.5 Hz probe field 
with amplitude of 0.02 Oe. Both exhibit pronounced non- 
linear responses, but with different forms most likely due to 
random field effects [l3] appearing for the more concentrated 
sample. 



germane, the samples respond very differently and the 
phase boundaries actually cross. 

Fig. [2] illustrates the pronounced sensitivity to fields 
applied parallel to the Ising axis and the evolution of 
the nonlinear response with x. xi^i) fo^" ^ = 0.167 can 
be described by a conventional power series expansion: 
X = Xi ~ 3x3^1 + 5x5 ^£ + • • '7 with all orders of the sus- 
ceptibility growing as the glass transition is approached 
from above. The longitudinal field dependence of the 
susceptibility for x = 0.198 also exhibits strong nonlin- 
earities, but has qualitatively different behavior. The 
parabolic xs at small h£ rolls over to a linear dependence 
at large field, consistent with a tendency towards the sin- 
gular linear behavior, attributed to random fields, seen 
for the disordered ferromagnet with x = 0.44 in the clas- 
sical low Ht, high T regime [17]. The dramatic change in 
X3, paired with the phase boundary crossing in Fig. [TJl, 
represents the major new discovery of our present work, 
and points to significant changes in the underlying physi- 
cal mechanisms arising from a small change in x near the 
onset of ferromagnetic long-range order. 

We show in Fig. [3] the temperature dependence of 
the linear and nonlinear terms in the susceptibility for 
X = 0.167 at Ht = 2.1 kOe. As expected, increas- 
ing orders of the susceptibility diverge increasingly more 
strongly refiecting the approach to a phase transi- 
tion. This is in accord with the results reported by Wu 
et al. [4] and in disagreement with the recent results of 
Jonsson et al. [23]. The discrepancy can be understood 
by looking at the different limits in which the system 
was examined. In Ref. |J] and in the present work, great 
care was taken to accumulate data in the hi ^ limit 
(±10 - 20 Oe about the peak). By contrast, the /iSQUID 
technique of Ref. [23] involved polarizing the system in 
a large longitudinal field (3 kOe, 15 times the scale of 
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FIG. 3: (color online) Linear and nonlinear susceptibilities for 
LiHoo.i67Yo.833F4 derived from data akin to Fig. [2^ and fit 
to X = Xi ~ ^X^hj + ^x^hj + • • •. In (b) we demonstrate the 
danger of not explicitly accounting for high order, non-linear 
effects. Fitting our data to a simple parabola out to 0.3 kOe 
or rapidly sweeping the field as in Ref. 0| (open circles) 
suppress the divergence of xs and mask the true physics. 



Fig. [2]), and then rapidly decreasing the field through 
zero at sweep rates of up to 50 Oe/s, corresponding to 
effective frequencies in this work of over 100 Hz. As can 
be seen from Fig. [H frequencies in this range are far from 
the equilibrium limit even in the presence of substantial 
transverse fields, suggesting that the results reported in 
(2^ do not capture the true physics. 

The differences in the underlying mechanisms that 
were observed in the nonlinear susceptibilities also 
emerge in dynamical measurements of the two concentra- 
tions. These can be determined by examining the spec- 
troscopic response of the system over several decades of 
/, as shown in Fig. [H for x = 0.198 at a series of Hf. This 
sample is slower than its x = 0.167 counterpart [3], where 
quantum fluctuations promoted by the off-diagonal ele- 
ments of the dipolar interaction more effectively speed 
the long-time relaxation. We characterize the approach 
to the spin glass from above by fitting the low frequency 
tail of Xi to a power law form, Tg{Ht) and Hc{T) are 
defined dynamically when a ^ and fluctuations occur 
on all (long) timescales. At this point xi grows loga- 
rithmically with /, where the onset frequency /o defines 
the fastest relaxation process available to the system and 
characterizes the quantum tunneling rate. 
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FIG. 4: (color online) Spectral response of the magnetic sus- 
ceptibility of LiHoo.i98Yo.802F4 at T = 0.05 K for a series 
of transverse fields. The spin glass transition is marked by a 
flat low- frequency response in Xi , corresponding to a logarith- 
mic dependence of Xi with onset fo and, via the fluctuation- 
dissipation theorem, 1// noise in the magnetization. 



The quantities derived from the nonlinear measure- 
ments and the spectroscopic response are combined in 
Fig. [51 Magnetic glass transitions for both x = 0.167 
and 0.198 are defined classically (a-d) and quantum- 
mechanically (e-h) by a ^ (c,g). When it is possible 
to reach the / ^ limit at modest Ht^ then a sharp, dy- 
namical feature in Xi (a,e) and a divergence of xs (b,f) 
also serve to define Tg{Ht). Both concentrations show 
similar behavior in the classical limit, but differ substan- 
tially in the quantum limit, underscoring the crossover 
in the underlying physics. The onset frequency for re- 
laxation, /o, follows an Arrhenius law, e~^/^^-^, in the 
classical limit where thermal fluctuations dominate (d) 
and a WKB form in the quantum limit (h) [lH. As a 
function of T with Ht = 0, fo for the two concentrations 
are indistinguishable. By contrast, at base temperature 
the JYt-dependent fo curves have substantially different 
slopes for the two values of x, with faster relaxation at low 
Hf in the x = 0.198 sample with the suppressed quantum 
glass transition. The slope seen for x = 0.198 is similar to 
what was previously observed for the x = 0.44 ferromag- 
net [11], suggesting that the random- field effects which 
play an important role in the dynamics of the ferromag- 
net are also significant in the x = 0.198 glass. 

We have verified that in low transverse fields the di- 
lute dipolar-coupled magnet, LiRoxYi-x^A^ can display 
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FIG. 5: Dynamical and non-linear signatures of the spin glass 
transition in the classical {Ht = 0, a-d) and quantum (T 0, 
e-h) regimes. The contrasting behavior of the two concen- 
trations with transverse field reflect the competing effects of 
quantum entanglement at small x and random fields at large 
X (see text). x=0.167 data in e and f was previously published 
in [4]. Onset of non-zero a in g corresponds to He. Solid lines 
in d and h are Arrhenius and WKB fits, respectively. Dashed 
line is the slope for x = 0.44 (111] . 



the static and dynamic signatures of a conventional spin 
glass. These results are in agreement both with early 
theory [21] and experiments of fifteen years ago, and in 
disagreement with incorrect (as shown above) conclusions 
drawn from interesting recent experiments [23[ where the 
same material was subjected to very strong and rapid 
perturbations away from equilibrium. While the static 
signature of the spin glass transition - a diverging non- 
linear susceptibility - seems to disappear for high i^t, the 
dynamical signature - the appearance of a fiat x"{f ~^ 0) 
- persists and indeed becomes sharper. This suggests that 
internal random fields 0, B 12 1 notwithstanding, there 
is a distinct quantum glass state that can be entered 
via a first order transition for which xs would 

not diverge. The striking new discovery that we make 
here is that the (quantum) critical field Hc{T 0) for 
this state is a non-monotonic function of x, with a lower 
value for x = 0.198 than for both the x = 0.167 spin 
glass and the x = 0.44 ferromagnet. We suspect that for 
X = 0.198 the random field effects seen near the Curie 
point for the x = 0.44 sample are important, and sup- 



press the magnetic glass phase. On the other hand, the 
glass-like state for x = 0.167 is more robust because of 
the larger quantum entanglement derived from the rel- 
atively greater population of antiferromagnetically cou- 
pled spins; these entanglement effects are known to play 
a major role in the dynamics of the 0.045 spin liquid 
0. As X is lowered even more, we land in the anti- 
glass spin liquid phase. The quantum spin glass then 
acquires new meaning as a valence bond glass [25[ of a 
type where there are multiple ways of drawing the bonds 
to construct pairs, while the spin liquid for x = 0.045 is 
a non-degenerate liquid with a unique pattern of valence 
bonds. 
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